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Abstract 

Recent cosmological observations suggest the presence of small but nonzero cos- 
mological constant A^^g ~ (2 x 10"'^ eV)^. It is an intriguing possibility that such a 
small cosmological constant is supplied by the potential energy density of an ultralight 
axion- like field (called as quintessence axion). If this axion couples to the electroweak 
gauge fields, its potential may be generated by the electroweak instanton effects. We 
calculate the axion potential assuming the supersymmetric standard model and obtain 
a surprising result that the induced energy density of the quintessence axion field is 
very close to the value suggested from the observations. 



1 Introduction 



The problem of cosmological constant is one of the most difficult problems not only in theory 
of quantum gravity but also in particle physics Its observational limit is some 120 orders 
of magnitude smaller than the natural unit of gravity. On the other hand, one might expect 
that zero point energies of dynamical fields would produce a cosmological constant of the 
order of the Planck scale. Even if there was a supersymmetry (SUSY) one would get an 
induced cosmological constant of order mgusY whenever the SUSY is spontaneously broken, 
where msusY is the SUSY-breaking energy scale. If one takes msusY — 1 TeV it exceeds the 
observational upper limit by 60 orders of magnitude. 

A number of attempts have been made to solve this problem, but no satisfactory solution 
has been found, so far [|l], |^, |^, ^, H, |], |^, |^. However, we expect that some yet unknown 
mechanism leads to an exactly vanishing vacuum-energy density {i.e. zero cosmological 
constant), since it seems easier to obtain zero than to obtain a minuscule finite value. 



Contrary to the above expectation, however, there is accumulating evidence [TQ] that 
the universe is indeed accelerating now. That is, the recent astrophysical observations 
strongly indicate the presence of a small but nonvanishing cosmological constant P, |l^ 

AL ^ (2 X 10-3 eV)^ (1) 

Thus, the situation becomes more puzzling and serious: we must explain not only why the 
cosmological constant Acqs is so small but also why it has the observed magnitude. 

We consider that any mechanism leading to the exactly vanishing vacuum-energy density 
may operate only at the true minimum of potential, since otherwise the inflation in the early 
universe does not take place. This argument has led ||Tl], |12|, to assume that the observed 
cosmological constant is the potential energy density of a dynamical scalar field A which 
differs from the ground state value. Namely, the potential V{A) has a stable minimum at 
which V = 0, but the field A has not yet reached the minimum point since the potential V 
is very fiat. This requires that the effective mass of A is not more than the expansion 
rate Hq of the present universe: 

m^<i7o - 2 X 10-33 eV. (2) 

The most natural candidate for such an ultralight field ^ is a pseudo Nambu-Goldstone boson 
of a nonlinearly realized global symmetry, which acquires its potential through instanton 



1 



effects [|Tll.Q There has been noted that such a global symmetry is explicitly broken by 
quantum gravity effects flSl. However, we postulate throughout this paper that those effects 
are sufficiently suppressed, since the dynamics of quantum gravity is not well known yet. 

The purpose of this paper is to calculate the potential of the pseudo Nambu-Goldstone 
boson A (called as quintessence axion) by using instantons of the standard electroweak 811(2)^ 
gauge theory.0 Surprising enough is that the electroweak instanton effects induce a correct 
magnitude of the potential energy density {i.e. the effective cosmological constant) given 
in Eq. (|ID. We consider that this result is very encouraging for the quintessential axion 
hypothesis. 

We assume the SUSY standard model throughout this paper, since it is a fascinating 
candidate beyond the standard model. In the minimal SUSY standard model (MSSM), the 
electroweak SU(2)l gauge interaction becomes asymptotic non-free above the SUSY particle 
threshold. We consider that the MSSM is an effective theory below the Planck scale and 
use the cutoff at the reduced Planck scale Mpi ~ 2 x 10^^ GeV to avoid the unphysical 
divergence. The S\J{2)l gauge couphng at the high-energy scale Mpi is a bit larger than 
that at low energies around the electroweak scale, which is one of important ingredients that 
induce the axion potential of appropriate order of magnitude. 



2 The Quintessence Axion Model 



= / d^e — -—w'^w^ + h.c. (3) 



We introduce the hypothetical axion superfield which couples to the electroweak SU(2)i 
field strength (a = 1, ■ ■ ■ , 3; a = 1, 2) as 

ill 

327r2 F4 

other than the MSSM particles, where F4 is the decay constant of the axion. This axion field 
$^ may arise from spontaneous breakdown of some global U(l)x symmetry at high energy 
of order F4, say the reduced Planck scale Mpi, or perhaps its origin can be traced to the 
string theory |T^. The superpotential Eq. (j^) possesses a global U(l)x symmetry generated 



by the shift $_4 $^ + iF^5, which has an anomaly of the SU(2)i gauge interaction. The 
quintessence axion field A is identified with the imaginary part of the lowest component of 
the superfield 

^ In the previous paper |Q, we calculated the potential V{A) assuming a hypothetical SU(7V) gauge 
symmetry broken at the Planck scale. 

The contribution to the axion potential from electroweak instanton effects has been discussed in Refs. 
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Table 1: U(1)r charges. Q, f/, L, and if^ represent the MSSM chiral superfields. S 

is the gauge-couphng superfield which is related to the SU(2)i gauge coupling as in Eq. (P). 
and V'^ are the super-covariant derivatives. 



In addition to the above global U(l)x, however, there is an accidental U(1)b+l symmetry 
in the MSSM which also has an SU(2)i gauge anomaly. There is, thus, a linear combination 
of U(l)x and U(1)b+l where the anomaly is canceled out, and the axion field remains exactly 
massless due to the presence of such an anomaly-free global symmetry. Therefore, explicit 
breaking of [/(1)b+l is needed to obtain the nonzero axion potential. The required [/(1)b+l- 
breaking terms are supplied by nonrenormalizable operators in the Kahler potential and 
in the superpotential.0 In this paper, we consider all possible higher-dimensional operators 
suppressed by the reduced Planck scale, as long as they are allowed by symmetries. Thus, the 
induced axion potential strongly depends on the symmetries at the Planck scale. We impose 
an U(l)/{ symmetry, since it is the most natural symmetry that guarantees the smallness 
of vacuum-expectation value (VEV) of the superpotential {W) ~ Mp\F [|1^, where F is the 
SUSY-breaking F term.Q The U(l)/j charges for the MSSM fields are given in Table |l[ All the 
MSSM interactions are invariant under this U(1)r symmetry, except for the constant term in 
the superpotential The U(l)i{ symmetry restricts possible forms of the (5 + L) -breaking 
nonrenormalizable operators in the Kahler potential and in the superpotential. In the next 
section, we estimate the axion potential by introducing U(l)j? symmetric nonrenormalizable 
operators. The breaking of U(l)j? is encoded in the nonzero gaugino mass and the U(l)ii;- 
breaking A terms, mi/2 — A ~ 1 TeV, in the MSSM sector. 



^ We assume that the i?-parity is exact throughout this paper. 

^ In any reahstic SUSY standard models, the SUSY breaking must occur at some intermediate scale, 
^ 10^" GeV. Then, to obtain the vanishing cosmological constant at the true minimum of the potential, 
the VEV of the superpotential should also be at the intermediate scale, lyV) ~ Mp\F, which is naturally 
ensured by the approximate i?-symmetry. 

^ This symmetry is the unique i?-symmetry which is consistent with the SU(5)gut and the neutrino 
masses generated through the see-saw mechanism |19{. 
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3 Estimation of the Axion Potential 



The axion potential is generated by the S\J{2)l instantons, which is given by integrating 
the size p of the instantons. The effects of the SU(2)i instantons are screened at distances 
larger than the electroweak scale by nonzero VEV's for the Higgs doublets, and hence we 
only have to integrate the region {h)~^ ■ Thus, we can calculate the axion potential in 
a manifestly supersymmetric manner in which all SUSY-breaking effects are described by 
VEV's of various superfield spurious [^ .p] 

Let us suppose that the effective Lagrangian at energy scale is given by 



Yi 



M, 



pi 



Y. 



m: 



di-3 



+ h.c. 



pi 



+ 



^'^7 + ] + h.c. 



(4) 



where operators Oi are the general operators composed of chiral superfields $r and $J, 
and supercovariant derivatives and Pq,, while Oj the holomorphic operators composed 
solely of $r's; di and dj {di,dj > 3) are the mass dimension of the operators Oi and Oj, 
respectively. We here consider the wavefunction renormalization and coupling constants as 
superfield spurious. 



Y 



s 



(1 + CiQ'^ + 



1 . e 



(5) 
(6) 
(7) 
(8) 

(9) 



where r/j and are the coupling constants; m,,, Cj, -Dj, Aj and mi/2 are the soft SUSY- 
breaking masses of order msusY — 1 TeV. 

The effects of one instanton and one anti-instanton of the SU(2)i can be seen by integrat- 
ing out all the superfields with vanishing VEV's. It is summarized in the following effective 



^ This approximation is, in fact, well justified, since the dominant contribution to the axion potential 
comes from the Planck-size instantons as shown below. 
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Lagrangian pO| : 



irefr(<f r, 2^", t)^, Zr, Y„ /i, /^t, Y„ F/, e^+^'; p) + h.c. 



+ 



J d'ee'('^'''''^^)w,s{<l>r,fi,Yf,p) + h.c. 



(10) 



where $r represents chiral superfields which have nonvanishing VEV's. The forms of the 
and Wes are constrained by the symmetries present in Eq. (H). The axion potential is 
obtained by substituting the VEV's for superfields, ($r), and spurious in Eqs. into 
Eq. ([T0|) and integrating Grassmann coordinates 6 and 6. 

To have the full axion potential, we further have to integrate the instanton size p. However, 
since it requires {B + L)-violating nonrenormalizable interactions to obtain nonzero potential 
for the axion, the integral may be divergent at the ultra-violet limit, p — 0. If it is the case, 
we use the cutoff at p ~ MpJ^ in the integral, and then the axion potential is dominated by 
instantons of the Planck size. 

Let us first consider the axion potential coming from the effective superpotential. Since 
both d'^9 and exp(— 871^(5* + $_4/F4)) carry the ^{1)r charge of —2 in the present model, 
Wes should have the U(l)/j charge of 4. Then, the form of Wes is restricted as 

W^, = p- (p^ (HMff [ipf^), J^J^) ' (11) 

where / is a holomorphic function which is constrained by the symmetries in the tree-level 
Lagrangian. From this argument only, we can find that the superpotential contribution to 
the axion potential V, 

V = A\(l-cos(4-]). (12) 



F 



A 



is bounded as 



Mpi 

~ (7 X 10-^ eV)^ (13) 

where the factor msusY is supplied by integrating the Grassmann coordinate 6"^ and we have 
used 0:2 (Mpi) ~ 1/23 to give the numerical value. The energy density obtained in Eq. ([T3|) is 
too small to explain the cosmological constant in Eq. (|l|) suggested from the observations. 



5 



Incidentally, the above conclusion does not change even if we introduce some new physics 
at the intermediate scale, as long as it does not change the beta-function of the SU(2)i. Let 
us extend, for example, the MSSM introducing nonzero neutrino masses through the see-saw 
mechanism |T^. In the region p^M^^, the effective superpotential contains the Majorana 



masses Mr for the right-handed neutrinos and the Yukawa couplings of them to the left- 
handed neutrinos and the Higgs field. Thus, the function / in Eq. ( |ri|) has an additional 
argument (pM/j) in this region. On the other hand, in the region p ^ the neutrino mass 
is described by the effective operator O = LLHuHu/Mji, so that the argument of / is given 
by l/(pM/j) in the region p^ M]^^. Therefore, the introduction of new physics always brings 
parameters smaller than order one as arguments of the function /, and does not enhance the 
energy scale A_4 of the axion potential. 

We next discuss the axion potential generated from the effective Kahler potential. The 
U(l)ij charge of Kes should be 2, and we find that the dominant contribution to the axion 
potential may come from the effective Kahler potential of the following form: 

K,, = p- {p&) f (k, J\ J\- 3 , J{- 3 , e^-^^^) , (14) 

where / is some function constrained by the symmetries of the model. We see that the 
resulting axion potential is dominated by the Planck-size instantons (p ~ ^pi^) which give 

^ e-^C(r/„A„A*)miusYMpi 

^ C(r/„ A„ A*) (1 eV)^ (15) 

where C(?7i, \j, A*) is a numerical coefficient depending on the coupling constants, rji, \j and 
A*, in the tree- level Lagrangian. Here, we have assumed the hidden sector SUSY-breaking 
scenario where the soft SUSY-breaking terms are nonvanishing at the Planck scale pT|.p] 
An explicit example of K^s is provided by considering the dimension-six and dimension-five 
operators O = QQU'^E'^ and O = QQQL in the Kahler potential and the superpotential, 
respectively. Then, together with the MSSM Yukawa couplings the Kes is written as 



= p--'{p&){ZaZE)-'{YQUHYLEH, 

~ 2 

(YoqeieA (Yqqql\ r f s+s^-j2rT^^^^A rifii 



^ The SUSY breaking induces the mass of the order of the SUSY-breaking scale mgusY — 1 TeV for SUSY 
partners of the axion field A. 
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Table 2: 11(1)^ charges. 



and we obtain the axion potential in Eq. ([15|). The potential for the axion field is given 
by summing up the contributions from all possible operators, but its order of magnitude is 
determined by the dominant contributions in Eq. (plsl). Note that the U{1)r symmetry is 
crucial to obtain the above result. If we do not impose it, we obtain a larger axion potential 
by a factor of Afpi/msusY- 

The numerical value given in Eq. (|1^) seems somewhat large to explain the cosmological 
constant in Eq. (|l|). However, it strongly depends on the size of various coupling constants 
in the coefficient C. In the next section, we argue that the coefficient C is much smaller than 
order one in a realistic SUSY standard model and find that in such a model the resulting 
cosmological constant is indeed very close to the value suggested from the observations. 



4 The Axion Potential in a Realistic SUSY Standard 
Model 

In the explicit example in Eq. (0), we have used the dimension-five operators O = QQQL 
to provide the required explicit breaking of U(1)b+l- However, this operator causes danger- 
ous d = 5 proton decays |Q. Thus, the coefficient of this operator should be sufficiently 
suppressed to evade too fast proton decay. The suppression of desirable amount is naturally 
attained by imposing the U(l)i? flavor symmetry to the quark and lepton multiplets. The 
11(1)^7 symmetry also provides realistic quark and lepton mass matrices through the small 



breaking parameter (0) /Mpi = e [^. Therefore, we here impose the U(l)i? symmetry and 
evaluate the size of C in the present model. 

We assign the U(l)i? charge consistent with the standard SU(5)gut, for simplicity. The 
charge assignment is given in Table 0, where lOj D {Qi, tJi, Ei} and 5^ D {Z)j, L^} are the 
quark and lepton chiral superflelds embedded in 10 and 5* representations of the SU(5)gut, 
respectively, and i = 1, ■ ■ ■ , 3 denotes the generation index. This charge assignment is moti- 
vated |2^ to explain the observed quark and lepton masses and mixings including the large 



mixing angle for the atmospheric neutrino oscillation [^. The U{1)f breaking parameter e 
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is determined as e ~ 1/17 to reproduce the quark and lepton mass matrices ||2^. It should be 
noted here that the present model is not necessarily consistent with the standard SU(5)gut, 
since the axion is supposed not to couple to QCD gluon fields. Therefore, we use the 
SU(5)gut representations only for a classification of the quark and lepton multiplets. 

The above 11(1)^7 has an 811(2)^ gauge anomaly, so that the factor exp(— 87r^(S' + $^/F4)) 
carries the U(l)i? charge of 10. As a result, the coefficient C{rii, Xj, A*) should be suppressed 
at least by the tenth powers of e to match the U{1)f quantum numbers. In the explicit 
example in Eq. (|16]), we find that C is suppressed exactly by the tenth powers of e if we use 
the operators O = Q1Q2UIEI, O = Q1Q1Q3L1, Q2Q2Q3L2, QsU^Hu and L^E^Hd to close all 
the fermion zero modes (see Fig. |I|). Thus, the dominant contribution to the axion potential 
is given by 

iv_4 _ e p ce mgugYiWpi 

^ c(^^j (IxlO-^eVr, (17) 

where c is an order-one constant, assuming that the magnitudes of all coupling constants are 
determined solely by the U(1)f charges. The resulting value of A^ is very close to the value 
Eq. (|l]) suggested from cosmological observations. This is very encouraging for identifying 
the axion field A with a quintessence axion. 

The axion potential Eq. (O) corresponds to the following axion mass: 



Pa 

10 . . .^18 ^ 2 



2 X 10^« GeV 



(7 X 10"^^ eV) 



Vl/17; V Fa 

This value satisfies the slow-roll condition Eq. (0) for the axion field as long as ^ 10^^ GeV.0 
If we set F4 ~ Mpi, this quintessence axion has an extremely small mass of order 10~^^ eV. 

We here comment on the numerical constant c. The constant c contains several am- 
biguities, but it is not far from of order unity. One possible source of ambiguities comes 
from loop factors, instanton measure and integration of the collective coordinate in instan- 
ton calculus, which may make c small. We give the counting rule for the tt^ factors in the 
Appendix. We find that these factors do not give extremely strong suppression of c; it is. 



* The condition Fj, ^ 10"'^* GeV does not necessarily need to hold, if the axion field lies precisely at the 
maximal point of the potential, 1 — cos((^) /FjC) ~ 2. 
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at most, only (I/tt^)"^ suppression. Since the axion potential receives contributions from all 
possible instanton diagrams, the number of diagrams may partially compensate the above 
suppression factors. Furthermore, there are ambiguities coming from the actual values of 
coupling constants and cutoff scale of instanton-size integrations. These ambiguities are also 
of a few orders of magnitude, so that it is reasonable to consider that the numerical constant 
c in Eqs. (0, [l^) is of order one. 



5 Discussion and Conclusions 

So far, we have calculated the axion potential assuming the particle content of the MSSM up 
to the Planck scale. However, if there exist some superheavy particles below the Planck scale 
Mpi carrying nontrivial quantum numbers of the S\J{2)l gauge symmetry, the SU(2)i gauge 
coupling (Mpi) at the Planck scale becomes larger and we may obtain larger potential 
energy of the axion. We now show that it is not the case if the masses for these superheavy 
particles are also determined by their U{1)f charges, and our conclusion on the induced 
cosmological constant remains unchanged. 

Let us introduce a pair of chiral superfields and ^ which belong to the vector-like 
representations of the S'U{2)l and have U(l)i;' charges q and q {q,q > 0), respectively. Then, 
these particles have a superheavy mass of order M^u — e^^^Mpi through the superpotential 
W = Msk"^^- The introduction of the superheavy particles changes the S'U{2)l gauge 
coupling at the Planck scale as 

a2(Mpi) «2(Mpi) 271 Msh' ^ ^ 

where Tsh denotes the Dynkin index of the SU{2)l representation to which belongs. This 
change enhances the instanton factor appearing in the axion potential as 

e -2(Mpi)_,g -2(Mpi) — PL I _ (20) 



VM 



SH 



On the other hand, the introduction of the superheavy particles also leads to additional 
fermion zero modes around the instanton configuration. To close these additional zero modes, 
we have to use appropriate operators which contain and \E' fields. As a result, numerical 
factors C{r]i, Xj, A*) appearing in the axion potential has an additional suppression factor of 
order e2(ij+g)TsH^ compensating exactly the enhancement in Eq. (^). The similar argument 
can also be applied for the case where the superheavy particle has a Majorana-type mass. 
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Thus, we find that the existence of the superheavy particles does not affect the order of 
magnitude of the axion potential, A^. 

It is interesting to note that the the gauge coupling unification scale can be raised up to 
the reduced Planck scale if there exist the superheavy particles and \&2 at the intermediate 
scale which belong to the adjoint representations of the SU(3)c and the SU(2)i, respectively 
[p^] . Even then, the above argument shows that the observed cosmological constant Eq. (|l]) 
is explained by the energy density of the quintessence axion. A, discussed in this paper. 
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A Counting tt^ Factors 

One might think at first glance that instanton amplitude is very much suppressed by loop 
factors; the diagram given in Fig. ^ for example, has many "loops" . However, the meaning 
of the "loop" is not trivial in instanton calculus, and the naive guess above is not true. In 
this appendix, we give a counting rule for vr^ factors appearing in instanton diagrams. The 
result Eq. (p6D shows that the diagram Fig. |l], for example, has only (l/vr^)^ suppression. 

We first count the vr^ factors analogous to the loop factors in ordinary perturbation theory. 
The instanton calculus are usually formulated in real Euclidean space-time, where the factors 
(l/47r^) and vr^ are provided by the propagator of nonzero modes and the angular integration 
of interaction vertex, respectively. Thus, we obtain the tt^ factors, 

where / and V denote the numbers of internal lines of nonzero modes and interaction vertices 
not including the instanton center, respectively, while L and C are the numbers of loops 
surrounded only by nonzero modes and connected parts after eliminating the instanton center, 
respectively. 

Another source for the vr^ factor is the instanton measure. The instanton measure con- 
sists of collective coordinate measure, super-product of zero-mode wavefunction norms which 
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comes from the super- Jacobian of the path-integral functional measure, Pauli-Villars regular- 
ization factor associated with zero modes, and the weight factor e J ^ evaluated at the tree 
level Among these, the norms of the zero-mode wavefunctions and the regularization 



yield tt^ factors. The zero-mode wavefunctions are given by the symmetry transformation 
of the instanton configuration and the Higgs scalar configuration twisted by the instanton 
background which contain no vr^ factor. Thus, the zero-mode wavefunctions do not contain 
TT^ factor. However, since the norms of zero-mode wavefunctions are given by integrating the 
wavefunctions, 

norm = (^J (i'^x (zero- mode wave function)^^ oc (7r^)2. (22) 

they have vr^ factors. Since there are ATq gauge-boson zero modes, 2Tg gaugino zero modes 
and 2Tr matter-fermion zero modes around the instanton configuration, the super-product of 
zero-mode wavefunction norms yields (l/7r2)~(^'^G'-2TG-2Tii)/2^ Together with the {1 / ^/2^^)'^'^'^ 
coming from Pauli-Villars regularization, the total vr^ factor from the instanton measure is 
given by 

1 

(23) 



7r2 

The rest of the sources for the vr^ factors are integrations of the collective coordinates. 
Integration of Grassmann collective coordinates does not yield any vr^ factor. Integration of 
the zero modes of gauge-group orientation gives the volume of coset space SU(A^)/SU(A^ — 2) 
X U(l) for SU(A^) gauge group, which has (tt^)^^^. The remaining is only the instanton-size 
integration. 

Combining all factors, we obtain the following instanton amplitude: 

J ^^^0^;;^ (^) n(external fields) Ab^-^-'p'm\\ (24) 

where b = 3Tg — Tfi and 4 — D is the sum of the mass dimensions of external fields in which 
various coupling constants are included through the spurion fields in Eqs. (H-^. We find 
that the instanton-size integration is always infra-red convergent because of the exponential 
factor (Higgs effect) while not always ultra-violet convergent. For the case of — 6 < 0, the 
integral Eq. (^4]) is fully convergent and we obtain the tt^ factor as 

I ^L-C+Ta-{N-l)-\{D~b) 



for D - 6 < 0. (25) 
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This reproduces the known result derived in asymptotic-free gauge theories [^. On the 



contrary, if D — 6 is positive, the divergence at the small scale must be cut off at pcut, which 
yields the following vr^-counting rule: 



I .L-C+Tr-(N-1) 



D-b 



.PcutAfpl^ 



for D - 6 > 0. 



(26) 
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Figure 1: One anti-instanton diagram generating the axion potential, e Together 
with one instanton diagram which gives e*"^/'^-^, we obtain the axion potential V = A^(l — 
cos(^/F4)) (see Ref. 
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